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NOMENCLATURE 


Ao  * A1 


f (x,t) 

F (5) 
n 

grad 


see  (2.51)  and  (2.52) 
airfoil  chord 

lift  coefficient,  see  (2.10) 

lift  curve  slope  at  zero  angle  of  attack 
(/  radian  or  /degree) 

moment  coefficient  about  midchord,  see  (2.11) 
see  (2.50),  (2.53)  and  Appendix  B 
2-D  vector  divergence  operation 
airfoil  deflection  mode 
see  (2.48) 

2-D  vector  gradient  operation 


t , 1 


Kn(z) 

K(x) 


unit  vectors  along  x and  y axes 

2-2—  , reduced  frequency 

00 

modified  Bessel  function  of  order  n , Ref.  19 
kernel  function 


L(x)  , L(x,t)  lift  distribution,  see  (2.22) 
p local  pressure 


Joo  > ►'oo 


free  stream  pressure  and  density 
dimensionless  pressure,  see  (2.8) 

pressure  on  upper  and  lower  surface  of  airfoil 

absolute  value  of  quantity  q 
1 2 

7poouoo  dynamic  pressure 
U00C 

— — Reynolds  number  referred  to  full  chord 


dimensionless  time 


£+  C + 


- 


L 


c 


T (x) 
n ' 


Chebyshev  polynominal  of  the  first  kind,  Ref.  19 
fluid  velocity 


velocity  of  point  x„  on  airfoil  surface 

s 


u_ 


Un(x) 


V 

w(x,t) 

(x,y)  or  x 

-*■ 

x„ 


free  stream  velocity 

Chebyshev  polynominal  of  the  second  kind.  Ref.  19 
perturbation  velocity  field 
downwash  velocity 

Cartesian  coordinates,  see  figure  1 
point  on  airfoil  surface 


a 

y 

V 

a 

T 

n 


angle  of  attack 
.57722,  Euler's  constant 
kinematic  viscosity 
Re/1* 

airfoil  thickness  to  chord  ratio 


streamfunction,  see  (2.15) 
\l/2 


o 1 - 


2ik 


Special  symbol 

r-|2 


2 2 
3 3 

— j j , Laplace  operation 

3x^  3y^ 
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I . INTRODUCTION 


In  recent  years  there  has  been  considerable  interest  in 
the  effect  of  the  boundary  layer  in  unsteady  and  steay  aero- 
dynamic theory.  Refs.  1-13-  A definitive  oscillating  airfoil 
experiment  is  described  by  Satyanarayana  and  Davis  in  Ref.  14. 
The  status  of  work  in  this  area  has  recently  been  reviewed  by 
McCroskey,  Ref.  15 , and  also  by  Williams,  Ref.  8.  The  import- 
ance of  boundary  layer  effects  in  transonic  flow  and  turbo- 
machinery has  been  discussed  by  Tijdeman,  Ref.  16,  and  Gostelow, 
Ref.  17.  For  the  non-lifting  "panel  problem"  the  principal 
role  of  the  boundary  layer  is  fairly  well  understood.  The 
dominant  effect  is  an  amplitude  reduction  and  phase  shift  of  the 
inviscid  pressure  distribution  due  to  the  momentum  defect  in  the 
boundary  layer.  The  corrections  are  of  the  order  of  the  ratio 
of  local  boundary  thickness  (or  some  inverse  power  of  Reynolds 
number)  to  the  length  scale  of  the  surface  disturbance.  The 
only  significant  effect  of  viscosity  in  the  non-lifting  problem 
is  to  establish  the  basic  turbulent  velocity  profile.  Beyond 
this  the  problem  may  be  considered  inviscid  and  the  so-called 
inviscid  parallel  shear  flow  model  may  be  used  to  calculate 
the  dominant  effect  of  the  boundary  layer.  Theoretical  results, 
both  steady  and  unsteady,  correlate  very  well  with  experimental 
data . 

The  status  of  the  lifting  problem  is  not  so  clear.  First 
of  all,  it  is  difficult  to  argue  that  inviscid  boundary  layer 
effects  are  important  at  all  for  a lifting  airfoil.  If  cor- 
rections are  of  the  same  order  as  in  the  non-lifting  problem 
(a  reasonable  supposition)  then  the  thickness  of  the  boundary 
layer  at  the  trailing  edge  must  be  at  least  as  large  as  the 
geometric  thickness  to  yield  a significant  effect.  For  "well 
designed"  high  Reynolds  number  airfoils,  such  is  not  the  case. 

It  is  for  this  reason  that  Dowell  and  his  co-workers  (Refs.  9- 
13)  and  the  author  (Refs.  1 and  2)  has  concentrated  primarily 
on  the  calculation  of  boundary  layer  effects  on  control  surfaces 


where  the  boundary  layer  thickness  may  be  comparable  to  the 
control  surface  chord. 

However,  If  we  examine  steady  state  experimental  data  on 
lift  curve  slope  for  various  airfoil  sections  (see  Section  II 
of  the  report  and  Ref.  18),  there  is  a significant  discrepancy 
between  the  measured  and  theoretical  values.  The  effect  of 
geometric  thickness  complicates  the  situation  even  further  by 
increasing  the  lift  curve  slope.  However,  at  high  Reynolds 
number  (1  to  10  million)  is  typically  10  to  20$  less  than 

the  flat  plate  theoretical  vtlue  of  2ir  . It  is  difficult  to 
imagine  anything  other  than  the  boundary  layer  that  could 
reconcile  theory  and  experiment.  Further  experimental  evidence 
for  the  importance  of  the  boundary  layer  in  the  unsteady  problem 
may  be  found  in  the  recent  work  of  Satyanarayana  and  Davis, 

Ref.  1^.  For  an  oscillating  airfoil  they  show  that  the  phase 
and  the  amplitude  of  the  surface  pressure  distribution  near  the 
trailing  edge  departs  from  the  inviscid  result  (plus  "Kutta" 
condition)  for  reduced  frequency  around  0.8.  The  experimental 
evidence  supports  the  view  that  boundary  layer  effects  must  be 
important  for  the  complete  airfoil  while  we  are  faced  with  the 
perplexing  result  that  Inviscid  parallel  flow  models  cannot 
predict  the  magnitude  of  the  correction  needed  even  for  the 
steady  state  lifting  airfoil. 

The  answer  to  the  dilemma  can  be  found  by  a closer  examina- 
tion of  the  steady  Incompressible  problem.  Most  textbooks  that 
deal  with  the  theory  of  lift  point  out  that  the  "Kutta"  condition 
is  an  emperical  condition  that  accounts  for  the  action  of  viscosity 
near  the  sharp  trailing  edge.  The  "total  lift"  on  the  airfoil 
is  due  to  this  viscous  effect.  Indeed,  it  was  one  of  the  great 
triumphs  in  the  theory  of  lift  that  the  singular  effect  of 
viscosity  was  recognized.  If  the  total  lift  is  determined  by 
the  action  of  viscosity  in  the  sub-boundary  layer  then  we  can 
see  where  we  have  perhaps  been  wrong  In  searching  for  "inviscid 
boundary  layer  thickness  corrections."  The  first  order  of  busi- 
ness should  be  to  calculate  the  way  that  viscosity  establishes 
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the  global  circulatory  flowfield  and  the  Reynolds  number 
dependence  thereof. 


In  this  report,  we  carry  out  the  above  program  within 
the  framework  of  linearized  "\iscous"  thin  airfoil  theory, 
f Starting  with  viscous  boundary  condition,  the  "Kutta"  condi- 

tion is  derived  in  the  limit  of  high  Reynolds  number.  The 
correction  to  the  Kutta  condition  is  found  to  be  of  order 
(1/ln  Re)  , and  is  much  greater  than  inviscid  boundary  layer 
f thickness  corrections.  In  fact,  the  viscous  correction  com- 

petes with  and  often  dominates  the  effect  of  geometric  thick- 
ness. We  also  derive  the  unsteady  viscous  airfoil  equation  and 
discuss  the  singularities  of  the  kernel  in  the  light  of  recent 
t experimental  work  in  Ref.  14. 


II.  TWO-DIMENSIONAL  VISCOUS  STEADY 
INCOMPRESSIBLE  THIN  AIRFOIL  THEORY 


Formulation  of  the  Problem 


Consider  a two-dimensional  flat  plate  airfoil  at  angle 
of  attack  a in  an  incompressible  fluid  (see  Figure  1). 


u % p 

* 00  * Koo 


u + U 


-c/2 


Figure  1.  The  Lift  Problem 


The  steady  state  incompressible  viscous  equations  are: 


where 


div  u = 0 


u*-—  = - grad  h + vV2u 
_ -*■ 

3x 


h = (P  - P00)/Pc 


(2.1) 


(2.2) 


(2.3) 


The  viscous  boundary  condition  is 


Lim  u(x)  = u (x  ) 

S 5 

x - xs 


(2.4) 


where  x„  denotes  the  position  of  the  airfoil  surface  and  ue 
its  velocity.  Disturbances  must  also  decay  at  infinity.  The 
problem  is  to  calculate  the  pressure  distribution  on  the  airfoil, 


.ait,  ■ 


t 


0 


We  choose  the  semichord  c/2  as  the  reference  length 

and  u as  the  reference  velocity.  For  reasons  that  will  be 
00 

apparent  later,  we  choose  the  Reynolds  number  based  on  the 
quarter  chord  as  the  measure  of  viscous  effects;  i.e.. 


u°°c  Re 

a = = T 


(2.5) 


where  Re  is  the  usual  Reynolds  number  based  on  the  full  chord. 
Next  we  introduce  a perturbation  velocity  field. 


v = u - lu^ 

(2.6) 

f 

where  for  the  moment  |v|  need  not  be 
less  problem  for  v becomes: 

small.  The 

dimension- 

div  v = 0 

+ +.9v  = _ 3P  + 
71  ^ 9x 

(2.7) 

o 

P - p„ 

» U2 
Moo  00 

(2.8) 

% 

Lim  v(x)  = us(xs) 
X + xs 

- i 

(2.9) 

We  denote  the  pressure  on  the  upper  and  lower 

surface  by 

% 

superscripts  + and  - respectively, 

section  lift  and  moment  coefficients: 

and  define 

the  usual 
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(2.10) 


(2.11) 


where  qM  Is  the  dynamic  pressure.  Clockwise  moments 
are  positive  and  is,  by  definition,  the  moment  coeffici- 

ent about  the  center  of  the  airfoil. 


B.  Formulation  of  the  Linearized  Viscous  Problem 

The  above  formulation  is  exact,  although  the  resulting 
problem  is  theoretically  intractable.  If  we  invoke  the  usual 
"no-slip"  boundary  condition,  the  fluid  mechanics  of  the  boun- 
dary layer  and  wake  are  exceedingly  complex.  The  effect  on  the 
lift  of  the  retarded  flow  in  these  regions  has  been  the  prin- 
cipal aim  of  recent  investigations  (Refs.  1-13)  that  start  with 
the  so-called  inviscid  parallel  shear  flow  model.  An  essential 
deficiency  of  the  parallel  shear  flow  model  (with  finite  surface 
velocity)  is,  however,  that  one  must  always  make  some  statement 
about  the  nature  of  the  flow  at  the  trailing  edge  (usually  the 
Kutta  condition)  or  in  the  wake  to  fix  the  inviscid  eigensolu- 
tion  (see  Section  II)  of  the  problem.  If  the  eigensolution  is 
omitted  the  total  lift  is  exactly  zero  and  the  flow  at  the 
trailing  edge  is  singular.  No  matter  what  details  of  the  inviscid 
boundary  layer  and  wake  flow  are  included,  the  situation  Is  the 
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same . 


The  first  order  problem,  then,  is  not  to  calculate 
inviscid  corrections  due  to  retarded  boundary  layer  and  wake 
flow,  but  to  find  out  how  the  "Kutta"  or  other  condition  is 
established  via  the  action  of  viscosity.  Most  text  books 
that  deal  with  the  subject  of  lift,  point  out  that  "viscosity" 
is  ultimately  responsible  for  the  lift.  Yet  the  aerodynamicist 
seldom  includes  viscosity  explicitly  in  the  lift  calculation. 

The  purpose  of  the  following  discussion  is  to  show  precisely 
how  viscosity  is  responsible  for  the  "Kutta"  condition  and  the 
steady  state  lift  on  an  airfoil  and  how  it  leads  to  the  first- 
order  correction  due  to  finite  Reynolds  number. 

Corrections  due  to  inviscid  boundary  layer  effects  are 
known  (Refs.  1-13)  to  be  of  the  order  of  the  boundary  layer 
thickness.  To  argue  that  these  corrections  are  small  compared 
to  viscous  effects,  we  pose  the  following  boundary  conditions. 
Suppose  that  we  fix  a symmetric  airfoil  at  zero  angle  of  attack. 
Then  we  "remove"  the  boundary  layer  either  by  symmetric  con- 
trolled suction  or  by  moving  the  surfaces  of  the  airfoil  at 
the  free  stream  velocity  u^  . For  theoretical  purposes  the 
latter  approach  is  favored  while  the  first  approach  is  exper- 
imentally more  practical.  Next,  we  place  the  airfoil  at  a small 
angle  of  attack  a . The  linearized  surface  velocity ,( 2. 9) » 
becomes , 

us  - - Ju^a  (2.12) 


Also,  the  linearized  fluid  equations  become 


div  v = 0 


lv  „ _ 3P  + _1_  y2- 
3x  3^  20 


(2.14) 


We  proceed  one  step  further  and  introduce  a streamfunction 
such  that 


u « ii  v ■ _ li 

3y  ’ 3x 


(2.15) 


Then  our  linearized  boundary  value  problem  becomes 


3 / 3U)  . 1 „2 . 

5?  - 57  + 75  v * 


. L.  / a±  + ± v2* 

Sx  3x  2o  * 


with 


^y(x,y  = 0±)  = 0 


*x(x,y  = 0±)  = a > ~i  < x < i 


(2.16) 


We  note  that  the  streamfunction  satisfies  the  classical  Oseen 
equation  and  the  pressure  is  a harmonic  function;  i.e.. 


3 n2. 

57  7 * 


57v''* 


(2.17) 


Also,  observe  that  the  quantity  ^ V2 ip ) is  a 

harmonic  function  that  is  conjugate  to  P . That  is,  P and 
A satisfy  the  Cauchy- Riemann  equations  (2.16).  Thus,  the 
complex  pressure  P' = P + iA  is  an  analytic  function  of  the 
complex  variable  z = x + iy  . One  may  compare  these  obser- 
vations with  classical  potential  flow  where  the  potential  and 
streamfunction  are  conjugate  harmonic  functions.  We  do  not 
pursue  the  development  of  these  concepts  in  the  present  report. 


The  novel  feature  about  (2.16)  is  the  presence 
of  the  viscous  terms  and  the  associated  viscous  boundary  condi- 
tions. In  the  following  section  we  use  these  equations  in 
conjunction  with  the  techniques  of  classical  thin  airfoil  theory 
to  develop  an  integral  equation  for  the  lift  distribution  on 
the  airfoil. 

C.  The  Viscous  Airfoil  Equation 

It  follows  from  the  boundary  conditions  (2.16)  that  ^ 
must  be  a symmetric  function  of  y so  that  P must  be  anti- 
symmetric in  y . We  develop  a solution  of  the  problem  by 
introducing  a distribution  of  pressure  doublets  along  that 
part  of  the  real  axis  occupied  by  the  airfoil,  i.e.. 


(2.18) 


(2.19) 


It  is  easily  shown  via  standard  techniques  of  classical  thin 
airfoil  theory  that 

Lim  P(x,y)  = ± ttQ ( x ) -1  < x < 1 


i 


and  so 


where 


Q(x)  = “ 2tt 


L(x)  = P"(x)  - P+(x)  , 


(2.21) 


(2.22) 


the  lift  distribution,  is  the  principal  unknown  of  the  problem. 
Next,  we  calculate 


3P  = _ _3 if 

3x  3x  2tt  J 


LU)  ^ In  R d? 


1 

3y  2 it  J 


L(c)  ^ In  R dC 


_L  (11  _ J, 

3y  1 3x  2a  v 


(2.23) 


where  we  have  used  the  first  of  equations  (2.15).  The  second 
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of  equations  (2.16)  leads  to  a similar  result  with  replaced 

by  . Thus,  we  integrate  to  obtain 

d X 


11  _ 

3x  - 2a 


ht 


A m R « 


(2.24) 


The  solution  of  the  last  equation  can  be  represented 
in  terms  of  fundamental  solutions.  First,  we  note  that  the 
right-hand  side  of  (2.24)  is  a solution  of  Laplace’s  equation, 
so  that 


,1 

I L(5)  In  R dC 


(2.25) 


is  a convenient  particular  integral.  We  add  to  (2.25)  a 
distribution  of  fundamental  solutions  of  the  homogeneous 
equation 


(2.26) 


The  appropriate  fundamental  solution  that  is  symmetric  in 
y is 

ip  = eox  Kq[c(x2  + y2)1/2]  (2.27) 

where  K is  the  modified  Bessel  function  of  zero  order  (see 
o 

Ref.  19).  The  fundamental  solution  characterizes  the  physical 
process  of  vorticity  diffusion  and  convection  and  is  needed  to 
satisfy  the  viscous  boundary  condition.  The  natural  viscous 
parameter  a is  the  Reynolds  number  based  on  the  quarter 
chord  (see  p.ll).  Thus,  we  can  write  the  general  solution  of 
(2.24)  in  the  form: 

4>  = I Q(C)eo(x-C)  Ko(aR)  d£ 

-1 


+ 


L(C)  In  R d£ 


(2.28) 


The  strength  of  the  viscous  solution,  Q(£)  , is  determined 
by  the  viscous  boundary  condition  (see  (2.16));  i.e., 


Llm  ty 
y ->■  o± 


Lim 

y -*•  0± 


1 

I 


3K 


QU)ea(x~C)  gy2-  (oR)  d£ 


if 


L(C)  d5 


so  that 


and 


= + ttQ(x)  t = 0 


Q(x)  = 


L(x) 

2it 


= — f1 

2 IT  J 


* - fp  I L(C) 
-1 


1 


In  R + ea(x_^)Ko(aR)  j d^ 


(2.29) 

(2.30) 


(2.31) 


Finally,  we  apply  the  second  boundary  condition  (2.16)  to 
obtain  the  "viscous"  airfoil  equation; 


i 

-1 


L(C)  K(x-5)  d?  = a 


(2.32) 


where  the  kernel  is  given  by 


K(x)  = 


= Ij-  In  | x | + e°X  Ko(a|x|) 


= x + aeaX  (Ko(alxl)  ' K! C a | x | ) sgn  x) 


. (2.33) 


If  one  follows  the  foregoing  steps  with  viscosity  set  to  zero, 
the  counterpart  of  (2.32)  is  the  classical  airfoil  equation. 


where  the  integral  is  a Cauchy  principal  value. 


(2.3*0 


The  importance  of  viscosity  can  be  illuminated  by  com- 
parison of  the  Cauchy  kernel  in  (2. 3*0  with  the  kernel  (2.33). 
We  observe  that  for  o|x|  0 the  kernel  (2.33)  becomes 


* - a In  + 0[a1 2|x|ln  c|x|]  (2.35) 

Thus,  the  viscous  kernel  only  has  a weak  logarithmic  singu- 
larity, no  matter  how  large  the  Reynolds  number  may  be. 

Now  it  is  well  known  that  the  solution  of  (2.3*0  is  not 
unique  (see  Section  III),  while  the  same  equation  with  a 
logarithmic  kernel,  i.e., 

1 

± f L (O  In  (x  - 5)  dC  - o (2.36) 

-1 

has  a unique  solution.  (See  Appendix  A and  Ref.  20.)  Thus, 
we  argue  that  by  retaining  a finite  Reynolds  number  in  (2.33) 
that  we  should  be  able  to  derive  the  unique  lowest  order  lift 
distribution.  We  show  by  a moment  technique  in  Section  II  E 
that  this  is  indeed  the  case,  and  the  Kutta  condition  is  the 
correct  lowest  order  "uniqueness  criteria"  to  apply  to  (2.3*0! 

Before  we  consider  the  solution  of  (2.33)  we  point  out 
another  essential  difference  between  the  viscous  and  inviscid 
kernels.  In  the  inviscid  case  the  kernel  has  upstream/down- 
stream symmetry.  The  airfoil  does  not  know  whether  it  is 
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flying  forward  or  backwards. 
(2.33)  is  highly  directional, 
expand  (2.33)  for  large  o|x| 


On  the  other  hand,  the  kernel 
To  illustrate  the  point,  we 
to  obtain 


K(x) 


ax  -°° 


ax  -*■  +°°  (2.37) 


The  upstream  viscous  effect  decays  exponentially  with  in- 
creasing Reynolds  number  while  it  decays  algebraically  as 
l//a | x | in  the  downstream  direction.  We  shall  see  that  the 
asymmetry  of  the  viscous  kernel  leads  to  removal  of  the 
"trailing  edge"  singularity  of  the  inviscid  solution  in  low- 
est order  in  accordance  with  the  Kutta  condition. 

D,  The  Inviscid  Solution  - A Brief  Review  and  Discussion 

The  solution  of  the  inviscid  airfoil  equation  (2.3*0 
for  an  arbitrary  downwash  distribution  a(x)  can  be  expressed 
in  closed  form;  i.e., 


L(x) 


(2.38) 


where  A is  an  arbitrary  constant,  the  "eigensolution. " 

For  the  special  case  of  constant  angle  of  attack  it  can 
be  verified  by  direct  substitution  and  use  of  known  relations 
for  Chebyshev  polynomials  (Ref.  19,  p.  783)  that 


L(x) 


A 2ax 

A - x2  A - x2 


(2.39) 


20 


I 


' 


i 

i 


[ 


t 


a 


t 

' 


- 


The  total  lift  depends  only  on  the  eigensolutlon;  i.e.,  using 
(2.10)  we  get 


C 


L 


ttA 


(2.40) 


The  Kutta  condition,  stated  in  the  form  that  "the  pressure 
loading  on  the  trailing  edge  must  vanish,"  fixes  the  constant 


A = 2a 


(2.41) 


and  we  get 


so  that 


L(x)  = 2a 


= 2na 


(2.42) 


(2.43) 


The  moment  about  midchord  is  independent  of  A ; i.e.,  from 

(2.11) 

1 

f 

CM  - - 4 I L(x)  x dx 


'll 

* f jd-d*. 

J.  ^77 


-1 


IT 

2 a 


(2.44) 


The  principal  point  to  be  made  is  that  it  is  the  eigensolution 
and  "empcrical"  uniqueness  criteria  in  the  form  of  the  Kutta 
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condition  that  determines  the  total  lift. 

We  also  call  attention  to  the  sensitive  nature  of  the 
integration  involved  in  the  inviscid  eigensolution.  Consider 
the  integral 


(x  - C)/l  - V 


-1  < x < 1 


(2.45) 


For  x very  close  to  unity,  we  note  that  the  integrand  in 
(2.45)  is  positive  over  the  entire  chord,  except  for  the  last 
segment  of  the  trailing  edge.  Thus,  the  integral  over  a 
vanishingly  small  segment  of  the  trailing  edge  must  balance 
the  integral  over  the  remainder  of  the  airfoil.  Obviously, 
the  same  statement  can  be  made  about  points  near  the  leading 
edge.  With  the  nonsymmetric  viscous  kernel ,( 2 . 33),  it  is 
evident  that  this  delicate  balance  is  destroyed  near  the  edges 
and  the  eigensolution  cannot  satisfy  the  viscous  counterpart 
of  (2.45).  This  observation  further  illustrates  the  very 
singular  way  that  viscosity  establishes  the  lift  on  the  airfoil, 

E.  Derivation  of  the  Kutta  Condition 

The  foregoing  discussion  suggests  strongly  that  the 
classical  Kutta  condition  is  implied  by  the  viscous  thin 
airfoil  equation  (2.32).  While  a rigorous  theoretical  proof 
is  lacking  and  a detailed  numerical  solution  of  the  problem 
has  not  been  carried  out,  we  can  show  that  this  is  indeed  the 
case  by  a moment  technique.  We  adopt  a trial  solution  of 
(2.32)  of  the  form  of  the  inviscid  solution  (2.39). 


A,. 

L(x)  = — ^ 


a - x2  A - x2 


(2.46) 


» 


To  determine  AQ  and  A^  we  require  that  the  first  two 
moments  of  (2.32)  be  satisfied.  That  is 


L(£)  P0(£)  d?  = 2a 


_1  f 

2v  I 


L(£)  F1(0  d£  = 0 


(2.47) 


where 


-1 

1 


'„<«>  - f *n 

-1 


K(x  - £)  d : , n = 0,1 


(2.48) 


and  K(x)  is  given  by  (2.33).  Substitute  (2.46)  into  (2.47) 
to  obtain  the  following  pair  of  equations  for  AQ  and  A1  ; 


1 

o 

o 

o 

o 

o 

, 

f A ^ 

A0 

it 

1 

o 
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O 

rH 

rH 

O 

, A1j 

2a^ 


v 0 J 


(2.49) 


where 


XF  (O  £m 


C o T r 


nm 


i_  f vi 

2V1  /iT 


d£ 


n 

m 


0,1 

0,1 


(2.50) 


The  solution  of  (2.49)  is 


A1 


2a 


'01 


111]  C 

,C10  00 


(2.51) 
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2a 


(2.52) 


A0 


(C11/C10) 

C01  ■ (C11/C10)  coo 


In  Appendix  B we  evaluate  the  coefficients  Cnm  asymptotically 
for  large  Reynolds  number.  The  final  results  are  summarized 
below 


'00 


'01 


'10 


'11 


— — (In  64Re  + y) - — 

4/tto  4/tto 


= - 1 + 


4/i To 


(In  64Re  + y - 4)  + — 


4/tto” 


-1—  (In  64Re  + Y)  + -2—  - -== 

4/170  4/lTO  /l7  0 

— — — (In  64pe  + y - 4) — + — — 

4/i To  4/ttF  3/tto  (2.53) 


where  y - .57722  is  Euler's  constant  and  Re  is  the  Reynolds 
number  based  on  the  airfoil  chord.  Thus  we  obtain 


A0  ~ + 2a 


1 + ( ~ - 2m)  /an  6 4 Re  + y)  + 0 ' ln  Re 


/Re  / J 


(2.5^0 


A1  ~ - 2a 


1 + 0 


In  Re 
k /Re  , 


(2.55) 


For  infinite  Reynolds  number  we  recover  the  inviscid  solution 
that  we  obtain  with  the  Kutta  condition;  i.e., 


Lim  TAq  = 2a 


Re  -*■  00 


A^  = - 2a 


(2.56) 


We  note  that  AQ  and  the  total  lift  is  determined  by  the  ratio 
of  two  viscous  terms  and  ClQ  that  vanish  asymptotically 
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at  the  same  rate  0(ln  Re/'/Re")  for  large  Reynolds  number. 
The  reason  that  AQ  tends  to  + 2a  as  opposec  to  - 2a  is 
that  the  principal  effect  of  the  viscous  boundary  layer  is 
felt  downstream.  The  singular  and  subtle  viscous  origin  of 
the  lift  is  thus  resolved. 


F.  Correction  to  the  Lift  Curve  Slope 

The  factor  in  square  brackets  in  (2.5^)  is  the  multipli- 
cative Reynolds  number  correction  factor  that  should  be  ap- 
plied to  the  lift  curve  slope  for  the  flat  plate  airfoil; 
i . e . , 


2tt 


20/3  - 2tt  \ 
In  64Re  + y ] 


rad 


(2.57) 


The  significant  fact  about  this  result  is  that  the  calculated 
viscous  correction  factor  is  0(l/ln  Re)  , and  is  large  com- 
pared to  boundary  layer  thickness  corrections  that  we  expect 

-1/2 

to  vary  as  some  inverse  power  of  Reynolds  number  (e.g.,  Re 
or  Re-1^).  The  positive  sign  of  the  correction  term  in  (2.57) 
is  somewhat  disturbing  in  that  we  would  expect  viscosity  to 
reduce  the  lift  curve  slope.  We  can  argue  in  favor  of  a 
smaller  CT  if  we  correct  for  the  very  important  effect  of 
leading  edge  thickness.  For  a thin  Joukowski  airfoil  the 
theoretical  lift  curve  slope  is  (Ref.  21) 

CT  « 2tt  | 1 + — rj  /rad 

La  ^ 3/3  }/  (2.58) 

where  t is  the  thickness  to  chord  ratio.  Thus,  CT 
increases  with  increasing  geometric  thickness.  Now,  within 
the  framework  of  our  present  linearized  theory,  thickness 
(actually  leading  edge  curvature)  has  a very  important  effect. 

It  removes  the  square  root  singularity  in  the  lift  distribu- 
tion at  the  leading  edge.  When  we  evaluate  the  integrals 
(see  Appendix  B)  Cnm  given  by  (2.53)  we  notice  that  the 
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correction  terms  in  parenthesis  result  from  integration  over 
the  trailing  edge,  while  the  remaining  constant  terms  result 
primarily  from  integration  of  the  singular  lift  distribution 
over  the  leading  edge.  If  finite  thickness  effects  are 
included,  we  can  argue  that  the  leading  edge  terms  should  be 
omitted  in  (2.53)  or  at  least  reduced.  The  resulting  lift 
curve  slope  becomes 


%.  ° 2’  ("  + 3/3  Tj  \ 1 ' ln  + Y 


Geometric  Kutta  factor 

thickness 

factor 


(2.59) 


where  we  have  also  included  the  correction  due  to  geometric 
thickness.  Thus,  the  effects  of  thickness  and  viscosity  oppose 
each  other  and  for  typical  high  Reynolds  number  thin  airfoils 
the  two  effects  are  of  comparable  order  of  magnitude. 

7 

For  example,  at  a Reynolds  number  Re  = 10  the  viscous 
correction  (Kutta  factor)  in  (2.59)  is  approximately  0.8l 
giving  a zero  thickness  lift  curve  slope  of 


C ^ = 0.0888/deg.  / Zero 

a ( thickness 


(2.60) 


or  approximately  a 20%  reduction.  For  an  airfoil  thickness  of 
10%  the  value  Is  increased  to 


CL  = 0.098/deg. 
a 

The  corresponding  result  for  Re  = 10^  is 


CL  s o.095/deg. 
a 


(2.61) 


(2.62) 


only  a 3%  difference.  While  the  Reynolds  number  correction  is 
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large  the  variation  with  He  is  very  small  at  high  values 
of  Re  because  of  the  logarithmic  dependence. 

The  net  effect  of  Reynolds  number  and  thickness  is  to 
reduce  the  lift  curve  slope  from  the  classical  thin  airfoil 
result  CT  = .1097/deg.  (or  2ir/rad).  The  values  given  by 
(2.61)  and  (2.62)  are  typical  of  measured  data  for  most  NACA 
airfoils  at  high  Reynolds  number.  In  Fig.  2 we  compare 
theory  with  experimental  data  (Ref.  18)  on  a 0015  airfoil. 

Since  the  leading  edge  radius  of  the  four  digit  series  air- 
foils varies  as  the  square  of  the  thickness,  it  is  not  certain 
what  geometric  correction  should  be  applied.  We  have  used  the 
10%  thickness  results  (see  (2.6l)  and  (2.62))  for  the  purpose 
of  comparison.  It  is  certain  that  the  Reynolds  number  correc- 
tion is  in  the  correct  order  of  magnitude  and  has  the  correct 
trend  for  Reynolds  numbers  greater  than  10^  . The  universal 
nature  of  this  correction  factor  is  apparent  from  the  more 
extensive  results  given  in  Table  1.  The  measured  's 

agree  remarkably  well  with  our  nominal  theoretical  values 
for  a wide  range  of  airfoils,  both  cambered  and  uncambered. 

7 

We  even  note  the  slightly  larger  values  of  CL  at  Re  = 10 

a 

G.  Concluding  Remarks 

We  have  derived  a logarithmic  Reynolds  number  correction 
to  the  classical  Kutta  condition  and  have  further  shown  that 
it  is  comparable  in  magnitude  and  often  exceeds  the  geometric 
thickness  effect  for  high  Reynolds  number  thin  airfoils. 

While  the  basic  logarithmic  dependence  on  Reynolds  number  is 
believed  to  be  correct,  it  must  be  pointed  out  that  the  magni- 
tude of  the  correction  term  should  be  calculated  more  accur- 
ately. In  particular,  the  solution  of  the  viscous  thin  airfoil 
equation, ( 2 . 32),  should  be  calculated  in  detail  for  large  but 
finite  Reynolds  numbers.  Also,  the  leading  edge  singularity 
and  possibly  the  trailing  edge  should  be  corrected  for  nonlinear 
geometric  thickness  effects.  Finally,  we  remark  that  the  unsteady 


Figure  2.  Variation  of  lift  curve  slooe  with  Reynolds  number 


linearized  viscous  problem  should  be  solved.  A formulation 
and  discussion  of  this  problem  is  given  in  the  following 
section. 
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III.  EXTENSION  TO  UNSTEADY  INCOMPRESSIBLE  FLOW 
A.  Basic  Equations 

The  basic  equations  for  the  unsteady  airfoil  problem 
can  be  derived  from  the  time  dependent  fluid  equations  fol- 
lowing the  steps  given  in  the  previous  section.  To  avoid 
unnecessary  repetition  of  detail,  we  simply  extend  the  basic 
results  of  the  previous  section  given  by  equation  (2.16) 

We  choose  the  ratio  as  the  basic  unit  of  time.  Then, 

the  unsteady  counterpart  of  (2.15)  and  (2.16)  can  be  expressed 
in  the  form, 


9P  = 9_ 
ax  3y 


3P 

9y 


(3.1) 


^y(x,y  = 0±,t)  = 0 

^x(x,y  = 0±,t)  = - w(x,t)  -1  < x < 1 (3.2) 

where 

D_  _ 3_  3_ 

Dt  “ 3t  3x  * (3.3) 


w(x,t)  - H , -1  < X < 1 (3.4) 

and  f(x,t)  is  the  airfoil  deflection  mode. 

The  steps  leading  to  the  integrodif ferential  equation 
(2.24)  of  the  previous  section  are  also  the  same,  so  that 
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where  the  basic  unknown  lift  distribution  L(x,t)  is  now  a 

function  of  time.  The  boundary  conditions  (3.2)  are  to  be  A 

applied  to  the  solution  of  (3.5). 


Review  of  Inviscid  Theory 


There  is  an  important  fundamental  difference  between  the 
steady  and  unsteady  problems.  In  the  steady  state  the  vortices 
shed  from  the  upper  and  lower  surfaces  of  the  airfoil  at  the 
trailing  edge  are  equal  and  opposite  in  sign,  so  that  a com- 
plete cancellation  of  the  vortex  wake  results.  In  the  inviscid 
theory,  this  means  that  all  singularities  of  the  streamfunction 
are  concentrated  along  the  airfoil  chord.  In  the  unsteady 
problem,  there  is  an  imbalance  in  the  vortex  shedding  process 
at  the  trailing  edge  and  a time  dependent  vortex  wake  is  formed. 
The  strength  of  this  wake  is  a significant  factor  in  determining 
the  total  unsteady  lift  distribution  and  is  again  determined 
completely  in  the  inviscid  problem  by  criteria  that  we  specify 
at  the  trailing  edge  (usually  the  steady  state  Kutta  condition). 
To  show  how  the  wake  affects  the  kernel  of  the  unsteady  airfoil 
equation  we  briefly  review  the  classical  unsteady  inviscid 
problem  (see  Ref.  22). 
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Consider  (3.5)  with  o = “ and  assume  simple  harmonic 
motion,  i.e. 


q(x,t)  =(Re  q(x)e-ikt)  (3.6) 


is  the  usual  reduced  frequency  based  on  the  semi-chord . 


Then, 


ik* 


- h f L<«>  I?  ln  R « 


A 


(3.8) 


We  integrate  (3.8)  from  upstream  infinity  to  any  point  x 
to  get 

.1 

t tc\  xe.  f ^ik(x-s)  

dS 

J 


hj  L(5>  «•  f 


Mx,y)  = i-  | L(£)  d£*  I eik(x_s)  ds-  hr  In  R 


(3.9) 


where 


r2  = (s  - O2  + y2 


(3.19) 


Finally,  apply  the  inviscid  boundary  conditions,  (3.2),  to 
obtain  the  incompressible  form  of  Poisso's  integral  equation 


h{ 


4-  l LU)  K( x - O d£  = - w ( x ) 


(3.11) 
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where 


K(x)  = 


= 9_  r 

3x^ 


ik(x-s)  ds 


A 

= i + lkeik*l  si 


iks 


x lke  J 
.00 


ds 


(3.12) 


The  Cauchy  principal  value  integral  can  be  expressed  in  terms 
of  sine  and  cosine  integrals.  For  present  purposes  we  need 
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only  point  out  that  the  second  term  has  a logarithmic  singu- 
larity. Simply  integrate  (3.12)  by  parts  to  get 


x 

K(x)  = i + ik  In  | x | - k2eikx  j e_iks  In  |s|  ds 

— 00 


(3.13) 


The  last  two  terms  in  (3.13)  are  due  to  unsteady  vortex  shed- 
ding in  the  wake.  Far  downstream,  the  wake  simply  oscillates 

ikx 

as  e . The  absence  of  wake  decay  in  the  inviscid  formu- 
lation is  due  to  the  fact  that  all  wake  effects  are  collapsed 
into  an  infintesimally  thin  sheet  of  vorticity  across  which 
discontinuities  in  tangential  velocity  can  occur. 

The  presence  of  the  Cauchy  singularity  in  the  kernel 
means  that  the  unsteady  problem  does  not  have  a unique  solution. 
By  invoking  the  Kutta  condition  one  can  obtain  the  well  known 
closed  form  solution  of  the  incompressible  problem.  If  we 
invoke  the  condition  that  no  vortex  wake  shall  form,  we  obtain 
a different  solution  of  the  problem  (Refs.  1 and  22).  The 
first  solution  is  consistent  with  the  steady  state  results 
obtained  in  the  previous  section  and  is  approximately  correct 
for  low  reduced  frequencies.  On  the  other  hand,  at  higher 
reduced  frequencies  one  might  expect  the  vorticity  shed  from 
opposite  sides  of  the  airfoil  to  self  annihilate  within  a very 
short  distance  downstream  of  the  trailing  edge.  Also,  the 
mechanism  whereby  vorticity  is  shed  into  the  wake  is  probably 
not  consistent  with  the  classical  steady  state  Kutta  condition 
as  the  frequency  is  increased.  The  experiments  of  Satyanarayana 
and  Davis  (Ref.  l4  ) show  clearly  that  the  classical  solution 
breaks  down  as  the  reduced  frequency  increases  towards  unity. 

Again  we  expect  that  an  explanation  for  the  discrepancy 
is  to  be  found  in  the  solution  of  the  viscous  problem.  While 
we  have  not  carried  out  a detailed  numerical  solution,  the 
basic  formulation  of  the  viscous  problem  and  the  structure  of 
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the  resulting  kernel  all  reinforce  the  basic  idea  that  linear 
viscous  effects  are  responsible  for  the  discrepancy. 


C . The  Unsteady  Viscous  Airfoil  Equation 


To  develop  the  viscous  airfoil  equation  we  consider  the 
fundamental  solution  of 


If  - lkf 


9 , x 

si  ln  r * — 

p 


(3. 1*0 


2 2 
x + y 


that  corresponds  to  a pressure  dipole  at  the  origin. 
(3.9)  we  have  the  inviscid  solution. 


From 


4>0(x,y)  = 


> 

■I 


,ik(*-s)  ds  9 ln  R 


(3.15) 


2 2 
R = s + y 


from  which  we  calculate  the  vorticity 


_ to  = v2ip  = 2tt  [eikx  H(x)  6 (y)] 

O O d X 


(3.16) 


= 2 ir  [6  (x)  6(y)  + ikeikx  H(x)  6 (y)  3 

The  inviscid  vorticity  is  a point  vortex  at  the  origin  plus 
the  usual  oscillating  wake.  Let 


* U>0  + 'J'1 


in  (3.1*0.  Then 


|f  - 1W 


h v2*'  ~ ch  Celkx  H(x)  6 (y)]  (3,17) 


The  Green's  function  of  the  homogeneous  equation  is 
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G(x,y)  = 2.  eax  KQ(fir) 


ft  = a (1  - 2 ik/a) 


1/2 


(3.18) 


so  that 


= jjd£  dn  ea(x_°  Ko(^R)  elk5  H(C)  6 (n  ) . 

o 

35  e°<K-5)  K0<aH)  fj  elks  H(C> 

f eik5  |-  eo(x"C>  K (8R)  35 

/ 3x  O 


00 
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(3.19) 


and 


a 

I 


ip  = / eik?  d£  [In  R + ea^x“^  KQ(nR)] 


(3.20) 


The  unsteady  viscous  kernel  is  simply  the  x derivative  of 
ip  evaluated  on  the  airfoil;  i.e.. 


00 

K(x)  = k Je 

o 


ikC  dC  [In  | x-C  | + ea(x-°  KQ(n|x-C| ) ] 


■ [In  | x | + eOX  Kn(H|x|)] 


+ ik  [In  | x | + eax  K («|x| )] 


(3.21) 
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d£  [In  | x-C | + ea(x_5)  K (n|x-?|)] 
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where  we  have  integrated  by  parts  to  obtain  the  final  result. 

The  terms  are  grouped  in  order  of  decreasing  singularity  of 
the  inviscid  kernel  (compare  with  (3.13)). 

We  first  note  that  the  Cauchy  singularity  is  removed  by 
viscosity  for  any  Reynolds  number  or  frequency.  Thus,  we 
expect  that  the  resulting  logarithmic  singularity  will  dictate 
a unique  solution  as  it  did  in  the  steady  state  case.  The 
logarithmic  singularity  in  the  unsteady  part  of  the  kernel  is 
also  removed  by  viscosity.  The  relative  importance  of  the  in- 
phase  (real)  and  out-of-phase  (imaginary)  singularities  of  the 
kernel  should  have  a strong  effect  on  the  phase  of  the  complex 
lift  distribution.  Higher  order  frequency  and  Reynolds  number 
effects  will  result  from  the  last  term  in  (3.21).  The  effect 
of  viscosity  on  the  unsteady  solution  is  two  fold.  First,  it 
is  responsible  for  the  basic  circulation  lift  with  an  import- 
ant Reynolds  number  correction  to  the  steady  state  Kutta  con- 
dition. Second,  it  affects  the  phase  with  which  vorticity  is 
shed  into  the  wake.  In  recent  experimental  work  of  Satyanarayana 
and  Davis  (Ref.  14)  it  was  observed  how  the  phase  and  amplitude 
of  the  surface  pressure  departs  from  inviscid  calculation  with 
the  steady  state  Kutta  condition  applied.  For  a reduced  fre- 
quency around  0.8  they  note  a discrepancy  in  the  phase  and 
amplitude.  For  reduced  frequency  near  unity  there  is  further 
steepening  of  the  pressure  distribution  near  the  trailing  edge. 
All  of  the  effects  they  observed  are  linear,  thus  indicating 
that  direct  viscous  action  is  responsible  for  the  discrepancy. 

It  is  believed  that  the  basic  physics  of  these  unsteady  phenomena 
is  contained  in  the  viscous  airfoil  equation.  Detailed  numerical 
solutions  must  be  carried  out  and  the  results  compared  with  the 
experiments  of  Satyanarayana  and  Davis  (Ref.  l4). 
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IV.  CONCLUSIONS 


The  principal  conclusion  of  this  work  is  that  the 
classical  eraperical  "Kutta  condition"  is  correct  in  the  limit 
that  Reynolds  number  approaches  infinity  and  is  a direct  con- 
sequence of  the  viscous  boundary  conditions.  The  correction 
to  the  lift  curve  slope  due  to  finite  but  large  Reynolds 
number  is  of  order  (1/ln  Re)  and  is  comparable  to  and  op- 
poses the  effect  of  geometric  thickness.  The  lift  curve 
slope  is  reduced  by  approximately  20%  for  Reynolds  number 
between  1 and  10  million.  Comparison  with  experimental  data 
confirms  this  calculation.  The  quantitative  effect  of  vis- 
cosity on  the  circulation  lift  greatly  exceeds  the  effects  of 
boundary  layer  thickness  or  momentum  defect. 

The  unsteady  viscous  airfoil  equation  is  derived  and  the 
singularities  of  the  kernel  are  discussed.  The  effect  of  vis- 
cosity is  to  remove  the  Cauchy  singularity  in  the  inviscid  kernel 
and  replace  it  by  a logarithmic  singularity.  The  lowest  order 
frequency  dependent  terms  in  the  kernel  also  has  a logarithmic 
singularity  that  is  removed  by  viscous  action.  Thus,  we  expect 
to  see  Reynolds  number  phase  variations  in  the  surface  load 
distribution,  as  the  frequency  is  increased.  This  conjecture 
is  in  accordance  with  recent  experimental  data  of  Satyanarayana 
and  Davis  (Ref.  14). 
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V.  RECOMMENDATIONS 

With  the  formulation  of  the  viscous  thin  airfoil  problem, 
we  have  opended  up  an  entirely  new  and  important  area  for 
research.  While  we  have  derived  some  very  significant  steady 
state  results,  we  have  only  begun  to  examine  the  unsteady 
problems.  Tu  further  develop  this  subject  the  following  steps 
are  recommended. 

A.  Steady  State  Incompressible  Theory 

1.  Develop  numerical  or  analytic  solutions  of  the  viscous 
thin  airfoil  equation  for  large  Re  . 

2.  Include  geometric  thickness  effects  In  the  viscous 
solutions . 

3.  Compare  lift  curve  slope,  center  of  pressure  and 
pressure  distributions  with  experimental  data. 

B.  Unsteady  Incompressible  Theory 

1.  Evaluate  the  kernel  of  the  viscous  unsteady  airfoil 
equation  and  develop  numerical  solutions  for  large 
Re  and  reduced  frequencies  up  to  order  unity. 

2.  Correct  the  solution  where  necessary  for  geometric 
thickness  effect. 

3*  Compare  results  with  experimental  data  of  Satyanarayana 
and  Davis  (Ref.  14). 

C.  Unsteady  Compressible  Theory 

Carry  out  a program  similar  to  A and  B with  the  additional 
feature  of  calculating  the  acoustic  far  field.  Both  surface 
pressure  distribution  and  the  acoustic  field  can  be  compared 
with  experimental  data,  e.g.,  those  of  Brooks  (Ref.  23). 


In  the  design  of  current  generation  aircraft,  it  is  often 

necessary  to  extrapolate  experimental  airfoil  data  in  the  range 

6 7 8 9 

10  < Re  <10  up  to  Reynolds  numbers  like  10  or  10  . 

Theoretical  results  of  the  type  that  would  be  obtained  with  the 

above  program  would  be  invaluable  for  establishing  intelligent 

extrapolation  procedures.  Reliable  theoretical  procedures 

would  also  preclude  the  necessity  for  building  costly  ultra- 

high  Reynolds  number  test  facillities. 
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APPENDIX  A 


A Uniqueness  Theorem 

Consider  the  homogeneous  integral  equation 

F(x)  = j g(£)  In  | x— C I d£  = 0 -1  < x < 1 

-1 


(A.  1) 


Is  there  an  cigensolution? 

Suppose  that  g(x)  Is  a solution  of  (A.l).  Then  the 
derivative  of  (A.l)  must  also  vanish;  I.e., 


0 -1  < x < 1 


(A. 2) 


Let 


g(x) 


n 


Tn<x) 


(A. 3) 


where  the  (x)  are  Chebyshev  polynomials  of  the  first  kind. 
They  form  a complete  set  of  function  on  the  interval  (-1  < x < 1) 
so  that  (A. 3)  is  a general  representation  of  the  solution.  Now 
substitute  (A. 3)  into  (A. 2)  to  get 

r1  t (o  dc 

) A_  i = - 7T  > A U , (x)  - o 

/ . n f > ■k  n n-1 

n^O  * A - C2(x  - O n=l  (A. 4) 

where  U^Cx)  the  Chebyshev  polynomial  of  the  second  kind. 

But  the  U^Cx)  are  a^so  complete  so  that  all  coefficients, 


A = 0 for  n-1 
n 
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(A. 5) 
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is  the  only  solution  of  (A. 2).  Is  it  a solution  of  (A.l)? 
Substitute  (A. 6)  into  (A.l)  and  choose  x = 0 so  that 


1 


The  integral  is  clearly  not  zero  so  that  Aq  = 0 and 

g(x)  = 0 -1  < x < 1 ( 


is  the  only  solution  of  (A.l) 
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APPENDIX  B 


Evaluation  of  Integrals  C 


Consider  the  integrals  of  Section  II,  (2.50); 


_ i 

2”J 


with 


— n _ o»i 

/i~~7?  m = 0,1 
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FnU)  = j xn  K(x  - C)  dx  n = 0,1 


(B.l) 


(B . 2 ) 


K(x)  = ^ [In  | x | + eox  Kq ( ct | x | ) ] 


(B.3) 


Evaluate  Fn  , n = 0,1 


Fo  = j Cln  lx"^l  + ea(x-C)  KQ(c | x-C | ) ] dx 


= m + ea(l-?)  K0[a(l-C)] 
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+ e“a(l+5)  K0[ct(1+C)] 

1 

-j  [In  | x-C | + ea(x_5)  K0(o|x-5|)]  dx 


(B.U) 
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The  coefficients  Cnm  can  be  expressed  in  the  form 


where 
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We  evaluate  each  of  the  integrals  for  large  a . Consider 


g(0  = j ea(x"C)  Ko(ct|  x-c| ) dx 


Let 


and  write 


t = cr(  x-O 


*0(1-0  , 

= ?J  K0(t)dt  + ?J 


g(0 


0(1-0 


e KQ(t)dt 


(B.13) 


(B. 14) 


For  a 00  the  first  term  is  uniformly  of  0(l/a)  or  smaller. 
The  principal  contribution  to  the  second  term  is  from  the  upper 
limit.  We  replace  KQ  by  its  asymptotic  form  to  get 


0(1-0 


g(C)  ~ oj  ^ + °(f)) 

0 

-Jp-  + o(i) 


dt  + 0 
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(B.15) 


Substitute  the  last  result  into  (B.ll)  and  (B.12)  to  get 

1 r1 
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(B.16) 


(B.17) 


We  note  that  both  and  Ag  are  primarily  a result  of 

integration  over  the  "singular"  leading  edge  lift  distribution. 
The  magnitude  of  both  terms  will  be  diminished  by  the  effect 
of  leading  edge  curvature. 

Next  we  consider  the  evaluation  of  integrals  A^  through 
Ajj  . We  observe  that  the  integrands  of  A^  and  A^  are 
weighted  near  the  trailing  edge,  while  those  of  A2  and  A^ 
are  concentrated  near  the  leading  edge.  Also,  the  integrands 
of  A2  and  A^  decay  exponentially  away  from  the  leading  edge 


while  the  integrands 

of 

A1 

and 

A^  decay  algebraically. 

Thus,  we  expect  A^ 

and 

A3 

to 

be  the  dominant  contribution 

to  the  coefficients 

Cnm 

• 

To  evaluate  A^ 

or 

A3 

we 

split  the  interval  of  inte- 

gration  at  £ = 1 - 6 where  a6  is  sufficiently  large  that 

K (a6)  can  be  replaced  by  its  asymptotic  expansion.  For  example, 

° 1/2 
we  could  choose  6 = 1/a  . Consider 
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The  first  integral  is  elementary  and  the  second  integral  is 
evaluated  in  Ref.  23.  The  final  result  is 
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— — (In  64  Re  + y) 
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(B.19) 


4l 


where 


Y = 0.57722 


(B. 20) 


is  Euler's  constant  and 


= Mo  = 


(B. 21) 


is  the  Reynolds  number  based  on  the  airfoil  chord.  To  evaluate 
A^  we  follow  the  same  steps  as  above.  The  result  is 


A^  = — — (In  6M  Re  + y - M)  + of-) 
3 4/™  W 


(B. 22) 


where  the  4 is  due  to  the  additional  power  of  £ in  the  inte- 
grand of  A^  . 

The  evaluation  of  A ^ or  A^  is  simple  because  the  only 
contribution  is  from  the  immediate  vicinity  of  the  leading  edge, 


Thus  we  have 


-1+6 


* . 1 f .-»<!+«  K0[n(It;)]  d*  + 0(e-2.4, 

2 2"  I a Jnrx  ( ’ 


OO 

V=  [ e'T  K0<t)  £ 
2tt/2o  J ° vT 


+ 0(e"2a5) 


-i-  + 0(e"2a6) 


4/Jr a 


Similarly 


4/rrcT 


+ 0(e“2a(S) 


(B . 2 3 ) 


(B.24) 


We  emphasize  that  both  A2  and  are  a result  of  integration 

over  the  leading  edge.  With  the  integrals  A^  through  Ag 


we  evaluate  C 


The  results  are  summarized  below: 


Summary 


Cnn  = — — (In  64  Re  + y) 
4/i To 


4/fra 


-1  + — — (In  64  Re  + Y - *0  + — — 
4/tto  4/fra 


— — (In  64  Re  + y)  + — — 

4/ncT  4/ttcj  /no 


Ch,  = — — (In  64  Re  + y - 4)  — + — 

1 4/fra  4/fTa 


4/fTa  3 /fra^  (B.25) 


